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A new method for the experimental study of topological effects 

in the quark-gluon plasma 
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A new method is presented for the quantitative measurement of charge separation about the reac- 
tion plane. A correlation function is obtained whose shape is concave when there is a net separation 
of positive and negative charges. Correlations not specifically associated with charge, from flow, jets 
and momentum conservation, do not influence the shape or magnitude of the correlation function. 
Detailed simulations are used to demonstrate the effectiveness of the method for the quantitative 
measurement of charge separation. Such measurements are a pre-requisite to the investigation of 
topological charge effects in the QGP as derived from the "strong CP problem". 

PACS numbers: Valid PACS appear here 



INTRODUCTION 

Topological charge fluctuations play an important role 
in the structure of the QCD vacuum Jj|. It manifests in 
the breaking of chiral symmetry, as well as in the mass 
spectrum and other properties of hadrons. These fluctua- 
tions can also lead to the formation of metastable vacuum 
domains, especially in the vicinity of the de-confinement 
phase transition, in which fundamental symmetries (V 
and/or CP) are spontaneously broken [0 i.e. the so- 
called "strong CP problem" . Experimental evidence for 
such topological fluctuations have been largely indirect. 

Recently, it has been suggested that direct experi- 
mental signatures of topological fluctuations could result 
from quark gluon plasma (QGP) [quarks liberated from 
hadronic confinement] subjected to an intense (hadron- 
scale) external magnetic field, via the so called chiral 
magnetic effect (CME) [], Q. In brief, topological charge 
fluctuations in the QGP leads to an axial anomaly or 
local imbalance between left-handed and right-handed 
light quarks. In an intense magnetic field, these quarks 
move along the field to create a net electric current which 
results in a separation of positive and negative electric 
charges in the field direction. Evidence for the chiral 
magnetic effect has been found in recent numerical lattice 
QCD calculations ||. An axial anomaly can also result 
from an anomalous global symmetry current in the hy- 
drodynamic description of the QGP [Q. This results in a 
modification of the hydrodynamic current by a term pro- 
portional to the vorticity of the fluid, and manifests also 
as a separation of positive and negative electric charges 
perpendicular to the reaction plane. Hereafter, we term 
this as the chiral rotation effect (CRE). 



MEASURING TOPOLOGICAL EFFECTS 
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FIG. 1. Schematic illustration of the chiral magnetic effect. 
Colliding nuclei (depicted as circles) are moving in and out 
of the page respectively. The Magnetic field (B) and system 
orbital angular momentum (L) are perpendicular to the reac- 
tion plane (plane spanned by the impact parameter, b, and 
colliding nuclei direction). Note that the overlap zone is not 
necessarily aligned with B or L. 



43 coupled low viscosity QGP ]_ 1| but also the strongest 

44 magnetic fields [orthogonal to the reaction plane] attain- 

45 able in the laboratory ||l5| . Consequently, the chiral mag- 

46 netic effect or the chiral rotation effect is expected to lead 

47 to a charge asymmetry in the distribution of particles 

48 emitted about the reaction plane (see Fig. Q). Exper- 

49 imental studies of such an asymmetry could provide an 
so important avenue for investigating one of the most im- 
5i portant problems of strong interaction theory. 



41 Collisions between heavy nuclei at the Relativistic 53 The STAR collaboration has analyzed data from re- 

42 Heavy Ion Collider (RHIC), not only create a strongly 54 cent measurements of Au+Au and Cu+Cu collisions (at 
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55 \/s~nn = 200 GeV) in search of this charge asymmetry 

56 with respect to the reaction plane. To do this they con- 

57 structed a correlator which used the emission angles of 

58 like-sign (+ + or ) and opposite-sign (-1 — ) hadron 

59 pairs. The correlator is defined by the event average 



2*, 



(1) 



AQ ~ C<++> + C7(~) - 2C(+-) 



(2) 



was used to test for a charge separation [about the re- 
action plane] of the kind suggested by the CME and the 
CRE, after an appropriate correction for dispersion of the 
reaction plane. 

A charge separation has been reported by the STAR 
collaboration Q. However, its mechanistic origin is 
still under intense debate | !^|. One reason for this 
has been the observation that the correlator used in the 
STAR analysis may be sensitive to several well known 
"background" correlations such as elliptic flow, jets and 
momentum conservation Therefore, it is im- 

portant to develop and investigate new correlators which 
can overcome many, if not all, of these deficiencies. 

A full study of topological effects in the QGP and its 
implications for the "strong CV problem" , will undoubt- 
edly require further detailed measurements focused on 
accurate experimental quantification of the dependence 
of charge asymmetry on particle species, particle pt, 
collision-system deformation, event centrality and beam 
collision energy. Here, we present a new experimental 
correlator specifically designed to aid such investigations. 

Our technique involves a multi-particle charge- 
sensitive in-event correlator C C (A5), which is expressed 
as a ratio of two distributions; 



C C {AS) 



N(AS csep ) 
N(AS cmix ) 



(3) 



The numerator is a distribution over events of the event 



89 averaged quantity AS CS ep defined as 



AS csep — {S p + ) (S n ) 



where 
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(St) 



X>in(A^_) 



(4) 
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n and p are the numbers of negative and positive hadrons 
[respectively] emitted about the observed event plane 
f bp (m = n + p is the charge hadron multiplicity for 
an event) and Atp = <p — ^ ep where <j) is the azimuthal 
emission angle of the charged hadron. 



96 The distribution AS cm i X in the denominator in Eq. 3, 

97 is obtained by making event averages in a slightly diffcr- 

98 ent way; That is, Eq. 5 is used to evaluate the averages 

99 (Sp) and (£^) for p and n randomly chosen hadrons 

100 (irrespective of charge) i.e. 



60 where 4> a ,4>p denote the azimuthal emission angles of any 

61 pair of hadrons, and ^2 denotes the azimuthal orienta- 

62 tion of the estimated second order event plane. The dif- 

63 ference 



AS r . 



(s h P ) - (s h n ) 



(0) 



There are several important features of the new cor- 
relator C C (AS). First, it is constructed entirely from a 
real event; hence, it is pure in event class (centrality, 
vertex, etc). Second, it is rather insensitive to the back- 
ground correlations which influence reliable extraction of 
the magnitude of the charge-separation correlation (see 
discussion below). In what follows, we use detailed sim- 
ulations to demonstrate the expected trends, as well as 
the efficacy of C C (AS). 



no SIMULATION METHODOLOGY 

in The response of C C (AS) to a charge-separation signal 

112 was tested via a detailed set of simulations tuned to re- 

113 produce observed experimental features. The simulations 

114 included the following major steps for each event. 
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• The event plane was chosen at random from 27r. 
Charged particles were then emitted with an az- 
imuthal distribution with respect to this reaction 
plane as: 



N(Aip) (x (1 + 2v 2 cos Atp) 



2v4 cos(4Atp) 
f2aisin(Ay>), 



(7) 
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where the Fourier coefficients v 2 and U4 are the 
observed magnitudes of elliptic and hexadecapole 
flow, and a\ is the charge-separation signal of in- 
terest. The number and px distribution of parti- 
cles were tuned to match the experimentally ob- 
served distributions. The reaction plane was then 
dispersed according to the experimentally observed 
dispersion for the centrality selection under study. 

• Neutral decay particles (e.g. A and Kq) were emit- 
ted with respect to the reaction plane according to 
their observed flow patterns. The decay kinemat- 
ics of these resonances were followed so as to obtain 
the daughter particle directions and momenta. The 
relative abundance of the decay particles were con- 
strained by the requirement that the simulated and 
observed positive-negative charge pair correlations 
(obtained by the standard event mixing method) 
were in agreement. 

• Jet particles were emitted with respect to the jet 
axes in a manner which was consistent with the 
observed two-particle jet correlations. 
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FIG. 2. N(AS)cse P and N(AS) C mix distributions for simula- 
tions performed with al > for all events. 
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FIG. 3. N(AS)cse P and N(AS)mi& distributions for simula- 
tions performed with al < for all events. 
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• All emitted particles were passed through an ac- 
ceptance filter specifically designed to take account 
of the detector acceptance and consequently, repro- 
duce the measured inclusive distributions for posi- 
tive and negative hadrons respectively as a function 
of pt- 

• The simulated events were analyzed as if they 
were actual experimental events. The correlation 
C C (AS) was evaluated for the selected range of pr 
using the detected particles and the dispersed re- 
action planes using Eqs. || - [| 
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RESULTS FROM SIMULATIONS 
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152 Simulations were performed for a broad spectrum of 

153 scenarios. Here, we show a representative set of results 

154 which lends insight into the detailed nature of C C (AS), as 

155 well as its sensitivity to different sources of background 

156 correlations. 

is? The distributions for N(AS) csep (solid circles) and 
us N(AS) m ix (open circles) are compared in Figs. [| and 

159 H for a\ > and a\ < (respectively) for all events. 

160 Fig. | shows that for oi > the distribution for 

161 N(AS) csep is shifted to the right when compared to that 

162 for N(AS) m i X . Similarly Fig. || shows that for a\ < 0, 

163 the distribution for N(AS) CS ep is shifted to the left when 
i6s compared to that for N(AS) m ix- Figs. [| and || show 

167 the respective C C (AS) distributions which result from the 

168 ratio of the distributions shown in Figs. |^ and |^. They 
«9 indicate sizable deviations from a flat distribution with 

170 positive and negative slopes respectively. Note that a flat 

171 distribution would be indicative of no charge-separation. 

172 C C (AS) distributions are shown in Figs. ^ and [?] for sim- 
i?3 ulatcd events in which (i) 51% of the events were gener- 
174 ated with oi > and the other 49 % with a\ < 0, and 



FIG. 4. Distribution for C C (AS) obtained from the ratio of 
the distributions shown in Fig. H. 

175 (ii) 49% of the events with a\ > and 51% of the events 

176 with oi < 0. In both cases, an asymmetric concave dis- 

177 tribution is obtained but with opposite asymmetry. 

178 The distributions for N(AS) csep and N(AS) m i X ob- 

179 tained for a simulation in which 50% of the events were 
lso generated with a\ > and the other 50% with a\ < 
lsi arc shown in Fig. ^. This choice was made to mimic the 

182 effects of local parity violation implied by current mod- 

183 els of topological charge generation in the QGP. For this 

184 scenario, Fig. ^ indicates that although the two distribu- 

185 tions are strikingly similar, N(AS) csep is slightly broader 

186 than N(AS) m i X . This is made more transparent in Fig.^ 

187 by the symmetric concave shape obtained for C C (AS) 

188 from the ratio of these distributions. 

189 To investigate the influence of "background" correla- 

190 tions from flow, jets and resonance decays, several sim- 

191 ulations were performed with these correlations turned 
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FIG. 5. Cc(AS) correlation function obtained from the ratio 
of the distributions shown in Fig. ^. 



FIG. 8. Comparison of N(AS) CS e P and N(AS) m ix distribu- 
tions for simulations performed with al < in 50% of the 
events and al > in the other 50%. 



1.01 

1.005 

C C (AS) 

1 

0.995 
0.99 



-0.5 



i i 



AS 



0.5 



1.005 



C C (AS) 



0.995 



0.99 



-0.5 





AS 



0.5 



FIG. 7. Cc(AS) correlation function obtained with al > in 
49% of events and al < in 51% of events. 
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FIG. 6. Cc(AS) correlation function obtained with al > in 
51% of events and al < in 49% of events. 
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FIG. 9. Cc(AS) correlation function obtained from the ratio 
of the distributions shown in Fig. ^. 



192 on or off. Fig. [h] shows the correlation function for a 

193 simulation in which only flow correlations are turned on 

194 for all events, i.e. «2,4 0, a\ = and resonance de- 

195 cays are turned off. The flat distribution indicated by 

198 Fig. [To] shows that C C (AS) is insensitive to flow. In 
we contrast to Fig. |o|, Fig. [H] shows a convex shape for 

199 C C (AS) which results from the charge correlations associ- 

200 ated with resonance decays which tends to bring opposite 

201 charges closer together in azimuth than on average. For 

202 this correlation function, the simulation was performed 

203 with flow on, oi = and resonance decays on, for all 

204 events. Since these correlations have an opposite infhi- 

205 ence on the shape of C C (AS) [compared to that for the 

206 charge separation signal] , it is important to have the rela- 

207 tivc abundances of decay particles properly incorporated 

208 into the simulations. This is ensured by requiring the 
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FIG. 10. Cc(AS') correlation function obtained for simulated 
events with only flow correlations. 
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FIG. 12. Cc(AS') correlation function obtained for simulated 
events with jets on, flow on, no resonance decay and ai=0. 



1.01 

1.005 

C C (AS) 

1 

0.995 



1.02 



1.01 

C C (AS) 
1 



0.99 



. a1=0.03 
o a1=0.06 



o 



o o 

O n O 



0.99 L 



-0.5 





AS 



0.5 



0.98 L 



-0.5 





AS 



0.5 



FIG. 11. C C (AS) correlation function obtained for simulated 
events with flow on, resonance decay on and ai = 0. 



FIG. 13. C C (AS) correlation functions obtained for simulated 
events with ai=0.03 and ai=0.06. 



209 standard two particle opposite charge correlations from 
2i» the simulations to match those for experimental data. 

212 The correlation function shown in Fig. [l^ was obtained 

213 for simulated events in which flow is on, a± = 0, resonance 

214 decays are off, but jets are turned on for all events. It is 

215 very similar to the flat distribution seen in Fig. [l^ and 

216 confirms the absence of any significant background cor- 
21? relations to C C (AS) from jets. Because the same event 
2ia is used to construct both N(AS) CS ep and N(AS) m ix mo- 

219 mentum correlation effects are also not expected to play 

220 any significant role. 

221 For an actual experimental correlation signal, the value 

222 of ai would be obtained by matching simulation to the 

223 observed correlation. The sensitivity of C C (A5) to the 

224 the parameter used to specify the magnitude of the 

225 charge separation a\ is demonstrated in Fig. |l3|. It 

226 shows that C C (AS) is responsive even to a relatively small 



227 charge asymmetry. 

228 It is important to stress that the method presented 

229 here is very general. For this study it has been applied 

230 to the correlations between charges in an event. However, 

231 our methodology can be applied to the investigation of 

232 correlations within any observed particle property. 



233 SUMMARY 

234 In summary, we have presented a new method which 

235 allows for good quantitative measurement of charge sep- 

236 aration about the reaction plane. Our method involves 

237 the formulation of a novel correlation function C C (AS) 

238 whose shape is concave only when there is a non-zero 

239 charge separation signal. The strength of C C (AS) is re- 

240 latcd to the parameter a\ which can be linked to a parity 



6 



24i viol ating signal. It is noteworthy that experimentally a 



242 rea dtion plane obtained from second order flow (vi) has 

243 a ±7r ambiguity. Introduction of such an ambiguity in 

244 the simulations makes all correlations symmetric i.e. the 



245 sensillivily Lo global pailly violation would be lost since 



the sign of «i cannot be deiei'mincd. However Hie sensi- 
tivity of the correlation function to the magnitude of a\ 



[io; 
in 



248 remains. Correlations due to flow, jets and momentum 275 [12 

249 conservation, do not influence the shape nor the magni- 

250 tude of C C (A5). Therefore, C C (AS) measurements could 

251 provide an important framework for investigating one of 

252 the most important problems of strong interaction theory 



253 - this "strong UP problem" 
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